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i Abstract 

. We prove that, in the non-extreme Kerr-Newman black hole geometry, the Dirac 

equation has no normalizable, time-periodic solutions. A key tool is Chandrasekhar's 
separation of the Dirac equation in this geometry. A similar non-existence theorem is 
^^O ' established in a more general class of stationary, axisymmetric metrics in which the 

Dirac equation is known to be separable. These results indicate that, in contrast with 
the classical situation of massive particle orbits, a quantum mechanical Dirac particle 
i must either disappear into the black hole or escape to infinity. 

O ■ 1 Introduction 

in 
o 



' It has recently been proved in |13| that the Dirac equation does not admit normalizable, 

0^ . time-periodic solutions in a non-extreme Reissner-Nordstrom black hole geometry. This 

Q I result shows that quantization and the introduction of spin cause a significant qualitative 

break-down of the classical situation, where it is well-known that there exist special choices 
of initial conditions for the motion of massive test particles which give rise to closed orbits 
. Indeed, the above theorem implies that the quantum mechanical Dirac wave function 
in the gravitational and electromagnetic fields of a static, spherically symmetric black hole 
^ ■ describes a particle which must either disappear into the black hole or escape to infinity. 

It is quite natural to ask whether this result is stable under perturbations of the 
background metric; i.e., if the non-existence theorem for normalizable periodic solutions 
of the Dirac equation remains true if the background metric and electromagnetic field are 
changed in such a way that the spherical symmetry is destroyed. This is precisely the 
question that we address in this paper. 

We are guided in our choice of a more general background geometry by the uniqueness 
theorems of Carter, Israel, and Robinson |^, |5|, from which we know that the most general 
charged black hole equilibrium state is given by the Kerr-Newman solution of the Einstein- 
Maxwell equations. Thus, in order to investigate the non-existence of periodic solutions 
in the most general black hole geometry, we have to study the Dirac equation in the Kerr- 
Newman background. The fact that this study is even possible rests on the remarkable 
discovery made by Chandrasekhar that the Dirac equation is completely separable into 
ordinary differential equations in the Kerr-Newman background geometry, even though the 
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metric is stationary and axisymmetric ||^]. In order to state our result, let us first recall 
that the Kerr-Newman solution is characterized by three parameters, namely its mass M, 
angular momentum aM and electric charge Q. We prove that the non-existence theorem 
for normalizable solutions of the Dirac equation remains true in the case of a rotating 
Kerr-Newman black hole provided that the angular momentum per unit mass and the 
charge are sufficiently small relative to the total mass of the black hole. Such black holes 
have both a Cauchy and an event horizon, and are thus referred to as non-extreme black 
holes. Specifically, we prove: 

Theorem 1.1 In a Kerr-Newman black hole for which a'^ + Q'^ < M^, the Dirac equation 
has no normalizable, time-periodic solutions. 



In Theorem 3.1, we give a similar non-existence result for the most general stationary, 
axisymmetric metric in which the Dirac equation can be separated by Chandrasekhar's 
method. The hypotheses of this theorem are stated as much as possible in geometric 
terms. While the solutions one obtains by imposing the Einstein-Maxwell equations on 
this metric are of limited physical interest, this theorem indicates that our non-existence 
result applies in a broader context (e.g. for a more general energy-momentum tensor). 
This paper is organized as follows. In Section E]^, we derive the separability of the 



Dirac equation in the Kerr-Newman geometry in a form slightly different from the one 



used by Chandrasekhar so as to recover the equations established in ||13| in the spheri 



cally symmetric limit, by letting the angular momentum parameter a tend to zero. In 



Section 2^, we work out matching conditions for the spinor field across the horizons. This 
gives rise to a weak solution of the Dirac equation in the physical region of the maximal an- 
alytic extension of the Kerr-Newman solution, which is valid across the Cauchy and event 
horizons. We then proceed in Section to establish the non-existence of time-periodic 
solutions. Just as in the spherically symmetric case, the crucial step consists in exploiting 
the conservation and positivity of the Dirac current to show that because of the matching 
conditions, the only way in which a time-periodic solution of the Dirac equation can be 
normalizable is that each term in the Fourier expansion of the spinor field in time and the 
angular variable around the axis of symmetry, be identically zero. While the regularity of 
the angular dependence of the separable solutions is manifest in the spherically symmetric 
Reissner-Nordstrom case, this is not so in the axisymmetric case treated in this paper. 
The regularity is therefore established in the Appendix. Finally, we consider in Section |3| 



the extension of Theorem 1.1 to more general stationary axisymmetric metrics. 



2 The Kerr-Newman Black Hole 

2.1 The Dirac Equation in Boyer-Lindquist Coordinates 

Recall that, in Boyer-Lindquist coordinates {t,r,'d,ip), the Kerr-Newman metric takes the 
form H 



ds"^ = gjk dx^x^ 

= — lot — a ^ , _ , , „„ , 



^ {dt - a sin2 ^ d^f -u(^ + dA - (a dt - (r^ + a^) dy,f , (2.1) 



where 

U{r,'d) = r^-Fa^cos^i? , A(r) = - 2Mr 
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and the electromagnetic potential is of the form 



a sin^ d dip) 



(2.2) 



The metric is singular at the origin r = and at the zeros of the function A. We shall 
consider the so-called non-extreme case > + Q^. In this case, A has two distinct 
zeros 



M 



M2 



and 



M + x/M2 



The two radii tq and ri correspond to the Cauchy horizon and the event horizon for the 
non-extreme Kerr-Newman metric, respectively. 

We briefly recall some elementary facts about the Dirac operator in curved space-time. 
The Dirac operator G is a differential operator of first order, 



G = iG\x)^ + B{x) 



(2.3) 



where B and the Dirac matrices are (4 x 4)-matrices. The Dirac matrices are related 
to the Lorentzian metric via the anti-commutation relations 



g'Hx) 1 



\{G\x), G\x)] 



]-(g^{x)G\x) + G''{x)G^{x) 



(2.4) 



The matrix B is determined by the spinor connection and the electromagnetic potential 
through minimal coupling. As such, it is determined by the Levi-Civita connection of the 
background Lorentzian metric (2.1) and the potenti al (|2.2| ). The Dirac matrices are not 
uniquely determined by the anti-commutation rules ( |2.4[) . The ambiguity in the choice of 
Dirac matrices adapted to a given metric is formulated naturally in terms of the spin and 
frame bundles fl^ . A convenient method for calculating the Dirac operator in this bundle 
formulation is provided by the Newman-Penrose formalism . More generally, it is shown 
in [12 1 that all choices of Dirac matrices satisfying ( |2.4D yield unitarily equivalent Dirac 



operators. Furthermore, in |12], explicit formulas for the matrix B in terms of the Dirac 
matrices G^ are given. In the following, we attempt to combine the advantages of these 
different approaches; namely, we first choose the Dirac matrices using a Newman-Penrose 



frame, and then construct the matrix B using the explicit formulas in [12 
We choose the so-called symmetric frame {l,n,m,m) of [Q, 



I 



n 



m 



m 



V2U\A\ 

V2ir\A\ 
1 



'2 2\ 9 

r^ + a^) 



d 



dt ~'~ dr 



d 



'2 2\ d ^ d 
-'^^^dt-'^d-r 



d 



d 



. . , d d 
la sm u — + -— + — 

2U V dt d-d sin-d d(f 

If B d id 

: — ^a smu — + -— — -— 

2U \ dt d-d sin-!? dip 



where e is the step function €{x) = 1 for x > and €{x) = —1 otherwise. (Because of its 
symmetry properties, this frame is somewhat more convenient than the Kinnersley frame 
used in [P|. Also, the notation with the step function allows us to give a unified form of the 
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frame both inside and outside the horizons.) The symmetric frame is a Newman- Penrose 
null frame; i.e. 

<l,n> = 1 , <m,rn> = —1 , 

and all other scalar products between the elements of the frame vanish. From this complex 
null frame, we can form a real frame (ua)a=o,...,3 by setting 

Mo = [l + n) , ui = -j={l-n) 

1 . — ^ 1 , _^ 

uo = (m + m) , U3 = (m — m) . 

This frame is orthonormal; i.e. 

9]kui4 = r]ab , V'"''uiu^b = 9^^ , (2-5) 

where r]ab = v'^^ — cliag(l, — 1, — 1, — 1) is the Minkowski metric. We choose the Dirac 
matrices 7°, a = 0, ... ,3 of Minkowski space in the Weyl representation 



^ I -1 ' ^ [ a 



where a are the usual Pauli matrices 

The 7" satisfy the anti-commutation relations 

ri"' = \ {7", 7"} • (2.7) 

We choose as Dirac matrices associated to the Kerr-Newman metric the following linear 
combinations of the 7*^: 

&{x) = uiix)^- (2.8) 



As an immediate consequence of (^^) and ( |2.7D , these Dirac matrices satisfy the anti- 
commutation relations ( |2.4| ). Next, we must calculate the matrix B in ( |2.3D . In [12|, it is 
shown that the Dirac matrices induce a spin connection D, which has the general form 

Dj = ^ - iE, - ieAj with (2.9) 

Ej = \ p{d,p) - ^ Ti{G^ V,G") G^Gn + TV(pGj V^G™) p , (2.10) 

where p = ^ejkimG'' G^G^G"^ is the pseudoscalar matrix, ejUm is the Levi-Civita symbol 
of curved space-time, and Aj is the electromagnetic potential. Using the spin connection, 
the Dirac operator (2^) can be written in the alternative form G = iG^Dj. Thus the 
matrix B is given hy B = G^{Ej + eAj). Since, in our context, the Dirac matrices G^ are 
linear combinations of the 7", the matrix p is simply the constant p = 7^ = i7'^7-'^7^7'^. 
As a consequence, the first and last summands in (2.10) vanish, and we obtain 

B = -^Tr(G™V,G")G^G™G„ + eG^Aj . (2.11) 
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Using Ricci's lemma, 

= 4Vj5"" = VjTr(G'" G") = TY((VjG™) G") + TV(G"^ (V^G")) , 
and the commutation rules of the Dirac matrices, we can simplify the trace term as 
Tr(G" VjG") G^GmGn 



Tr(G'" V,G") [51 G„ - 5'iGm + G^ gmn + iei^np YG^ 
-8VjG^ + ie^'-'^P TV(G™ V,G„) 7'Gp . 



(2.12) 



Since the Levi-Civita connection is torsion-free, we can replace the covariant derivative in 
the last summand in (p.l2|) by a partial derivative. We thus conclude that 



B 



{g denotes as usual the determinant of the Lorentzian metric). This formula for B is 
particularly convenient, because it only involves partial derivatives, so that it becomes 
unnecessary to compute the Christoffel symbols of the Levi-Civita connection. Next, we 
substitute the derived formulas for G^ and B into (|2.3| ) and obtain for the Dirac operator 



G 







e(A_)_a_ 
-P- 
1 /. d 



a+ /3+ 

/3_ e(A) a_ 
-^5+ 
a+ 

cot a sin 

+ i + 



with 



1 



U \ 

VuW\ 



d 



2 

cott? 



2U 
a sini? 
2U 



(r — ia cos i?)^ ± 
(r + ia cos ■!?) I it 



• ^ 1 

a sm w — + 



(2.13) 



d 



dt sint? 

.a Id 

a sm?? — + — 
at smxf o(p 



■ 1 1 2\ d . d ^ 
[lir +a) — +ia — + eQr 
\ at dip 



|A| 
IT 



d r-M 
or 



2A +^ir-iacos^) 



a± 



VuW\ 



[ilr +a) — +ia — + eQr 
V ot dip 



1^ 

IT 



d r — M i 



The four-component wave function ^' of a Dirac particle is a solution of the Dirac equation 



(G - m) ^ 



v 



— m 







P+ 





—m 


P- 


e(A) Q 


e(A) a_ 




—m 





-P- 







— m 



(2.14) 



For a ^ 0, the Kerr-Newman metric goes over to the spherically symmetric Reissner- 
Nordstrom metric. Before going on, we explain how one can recover the Dirac operator 
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of [13| in this limit: First of all, in |13] the Dirac (and not the Weyl) representation of 
the Dirac matrices is used. Furthermore, instead of working with orthonormal frames, 
the Dirac matrices are constructed in by multiplying the Dirac matrices of Minkowski 
space in polar coordinates 7*, 7*", 7'', 7*^ with appropriate scalar functions. Because of 
these differences, the operator C^^^ obtained from ( |2.13D in the limit a ^ coincides with 
the Dirac operators in only up to a unitary transformation. More precisely, 



we have, in the notation of |13], 

^in/out ^ ^a=0 




V{^,^) = —{ ^ ^ ) + + + exp(^-i|a^j exp(^-i-a 

We point out that the unitary transformation V is not continuous, as it changes sign along 
the line <p = 0. Since in the limit of flat Minkowski space, G°^^ goes over to the usual 
Dirac operator in the Dirac representation, it is natural to assume (as in [||]) that the 
Dirac wave function corresponding to cWout -g continuous (or even smooth away 

from the poles "i? = 0, vr). According to the above transformation of the Dirac operators, 
the wave function ^ in ( p.l4| ) is given by = y-l ^in/out ^ ^^^^ 

we conclude that it must 

satisfy the boundary conditions 

lim ^(t,r,'d,{p) = — lim ^{t,r,'d,(p) . 

As a consequence, we must for the eigenvalues of the angular operator id^ consider half 
odd integer eigenvalues (cf. (|2.17D below). 



It is a remarkable fact that the Dirac equation (2.14) can be completely separated 
into ordinary differential equations. This was first shown for the Kerr metric by Chan- 
drasekhar [0, and was later generalized to the Kerr-Newman background |l9[. We 
shall now briefly recall how this is done, since we will need the explicit form of the corre- 
sponding ordinary differential operators for the proof of Theorem |1 . 1|. We closely follow 



the procedure in Q. Let S{r,'d) and T{r,'d) be the diagonal matrice^ 

S = |A|4 diag (^(r — cosi?)2, (r — ia cost?)2 , (r + cosi?)2 , (r + M cost?)2^ 
F = — i diag {{r + ia cos —{r + ia cost)), —{r — ia cost}), {r — ia cos 

Then the transformed wave function 

4/ = S"^ (2.15) 

satisfies the Dirac equation 

rS{G-m)S-^^ = . (2.16) 



This transformation is useful because the differential operator ( p. 16 ) can be written as a 



sum of an operator TZ, which depends only on the radius r, and an operator A, which 
depends only on the angular variable i9. More precisely, an explicit calculation gives 

rS{G-m)S-^ = TZ + A 



^We mention for clarity that this transformation of the spinors differs from that in by the factor 
|A| 4 in the definition of S. Our transformation simplifies the radial Dirac equation; moreover, it will make 
the form of our matching conditions easier. 
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with 



/ imr \/W\T^+ \ 

-imr e(A) V]A[P_ 

e{A)./\A\V_ -imr 

^ y/\A\'D+ imr j 

^ —am \ 

am cosi) — £_ 

£+ —am CDS'!? ' 

\ — £_ am cos?? J 



where 



P± - — T - 
dr A 



- d coti? 



' . , d Id 
a sm-d — + - — - — 
at smv Oip 



Now for ^ we first employ the ansatz 



a; G IR, A; e Z 



(2.17) 



with a function which is composed of radial functions X-i-{r) and angular functions 
Y±{-d) in the form 

/ X_(r)y„(T9) \ 
X+(r) y_(T9) 

V x„(r)y+(^) y 

By substituting ( ^.171) and ( 2.18 ) into the transformed Dirac equation ( 2.16| ), we obtain 
the eigenvalue problems 



^{r,'d) 



(2.18) 



whereby the Dirac equation ( p. 16 ) decouples into the system of ODEs 



/ \/| A| imr — A 

I -imr -A e(A)V]A|P_ 

£+ —am cos ?? + A 
am cos 19 + A — £„ 



y_ 







where T>± and reduce to the radial and angular operators 



V — ^ ± ^ 
dr A 



1 



d 



coti9 



(r^ + a^) + I^A; + ^ ) a + eQr 
au! sint? + 



2, 



sini? 



(2.19) 

(2.20) 
(2.21) 

(2.22) 
(2.23) 



We mention that in the limit a ^ 0, Chandrasekhar's separation of variables cor- 
responds to the usual separation of the angular dependence in a spherically symmetric 
background. The angular operator A then has the explicit eigenvalues A = ib(j + ^) with 
j = ^, |, . . .. The regularity of the eigenfunctions of the angular operator A for general a 
is established in the Appendix. 
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2.2 Matching of the Spinors Across the Horizons 

Let us consider the Dirac wave function = S^^"^ with ^' according to the ansatz 
( |2.17| ), (|2.18 ). Then the Dirac equation ( 2.14| ) separately describes the wave function in 



the three regions r < tq, rg < r < ri, and r > ri; for clarity, we denote ^' in these 
three regions by "ifM, and ^I'o, respectively. Since the ODEs ( 2.20| ),( p^ ) and the 



transformation S are regular for r {0, ro,ri}, the functions ^m, and are 
smooth. However, the difficulty is that the coefficients in the Dirac equation have poles at 
r = vq and r = ri. As a consequence, the Dirac wave function will in general be singular 
for r — > rg^. Furthermore, is not clear how to treat the Dirac equation across the horizons. 
In this section, we shall derive matching conditions which relate the wave functions inside 
and outside each horizon. For the derivation, we shall first remove the singularities of 
the metric on the horizons by transforming to Kerr coordinates. In these coordinates, we 
can also arrange that the Dirac operator is regular. This will allow us to derive a weak 
solution of the Dirac equation valid across the Cauchy and event horizons. In the end, we 
will transform the derived conditions back to Boyer-Lindquist coordinates. 

First, we must choose coordinate systems where the metric becomes regular on the 
horizons. One possibility is to go over to the Kerr coordinates («+, r, given in 

infinitesimal form by Q 

-\- O?' CL 

duj^ = dt + — ^ — dr , dip+ = dip + — dr . 

Alternatively, we can choose the coordinates (n_, r, with 

7^2 _j_ q2 ^ 

du- = dt — dr , d(fi- = dip dr . 

A ^ ^ A 

The variables and u„ are the incoming and outgoing null coordinates, respectively. 
Along the lines u± = const, the variables t and r are related to each other by 

dt = =F : dr . 

A 

By integration, we see that, for r rg, we have t — > =pcxD. Similarly, lim,.^ri t = ^oo. 
The fact that t becomes infinite in the limit r rQ/i means that the Kerr coordinates 
describe extensions of the original Boyer-Lindquist space-time, whereby the Cauchy and 
event horizons have moved to points at infinity. More precisely, the Cauchy horizon 
corresponds to the points (r = rg, n_|_ = oo) and (r = rg, U- = — oo); the event horizon 
is at (r = ri, n+ = — oo) and (r = ri, u_ = cxd). In other words, the chart (n+, r, ??, 
extends the Boyer-Lindquist space-time across the points (r = rg, t = — oo) and (r = 
ri, t = oo), whereas the chart {u-,r, (/9_) gives an extension across (r = rg, t = oo) and 
(r = ri, t = — oo). 

We next work out the transformation of the wave functions from Boyer-Lindquist 
to Kerr coordinates. We first consider the chart (n+, r, -i?, (/94.). The transformation of 
the Dirac equation consists of a transformation of the space-time coordinates and of the 
spinors. For clarity, we perform these transformations in two separate steps. Changing 
only the space-time coordinates transforms the Dirac matrices to 

QU+ ^ Qt QT _h ^ Z ^ / U _ i\ 

dt dr VC/ \/UW\ 
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e(-A) 7°) 



1 



VuW\ 



(7° 



where is the Heaviside function Q{x) = 1 for x > and Q{x) = otherwise. The 
matrices G"+ and are singular on the horizons. Therefore we transform the spinors 
and Dirac matrices according to 



with 



V{r) 



V{r) ^ 






V{r) G^ V{r) 



i(|Ari + |A|i)l-l(|Ari 



lAli 



(2.24) 



(2.25) 



The transformed Dirac matrices are 



&+ = 


asin"!? 


7^ 


Vu 


& = 


1 






& = 


' 7^ 


G\ = 


1 




sin?? V 





+ a? 



U 



(7° - 7^) 



(1 - A) 7° + (1 + A) 7= 



7^ + 



U 



Now the Dirac matrices are regular except at the coordinate singularities = 0, vr and 
at the origin r = 0. The anti-commutation relations ( |2.4| ) allow us to check immediately 
that the metric is indeed regular across the horizons. The transformed Dirac operator G 
can be constructed from the Dirac matrices G^ with the explicit formulas (2.11) and ( |2.3| ) 
(these formulas are valid in the same way with an additional tilde, because the matrices 
G^ are again linear combinations of the Dirac matrices 7-' of Minkowski space). From this, 
we see that all the coefficients of the Dirac operator G are regular across the horizons. 



According to the transformation (2.24) of the wave functions, the Dirac operators G and 
G are related to each other by 

G = VGV~^ . 

Since the operator G is regular across the horizons, we can now study the Dirac equa- 
tion on the event and Cauchy horizons. We denote our original wave functions trans- 
formed to Kerr coordinates by ^ and ^o- They are smooth in the regions r < tq, 
rQ < r < ri, and r > ri and satisfy the Dirac equation there. However, they may have 
singularities at r = rg and r = ri. Let us assume that ^ := ^'/-|-^'jv./-|-^'o is a generalized 
solution of the Dirac equation across the horizons. In order to analyze the behavior of the 
wave function near the Cauchy horizon, we write ^' in a neighborhood of r = ro in the 
form 
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and substitute into the Dirac equation (G — m)^ = 0. Since ^ is a solution of the Dirac 
equation for r 7^ rg, , we only get a contribution from the derivative of the Heaviside 
function, i.e. in a formal calculation 

To give this distributional equation a precise meaning, we multiply the above formal 
identity by a test function r/(r) and integrate, 

0=/ r?(r)5(r-ro) (7° + 7^) (^Af(^i+,r,i9,v?+) - ^'/(n+,r,i9,v?+) . (2.26) 

By choosing a function r/(r) which goes to zero sufficiently fast for r — > tq, we can make 
sense of this integral, even if (7*^ + 7^)(^'m — ^ i) is singular in this limit. For example, 
we can choose r/ as r/ = (1 + 1(7*^ + 7'^)(^b — where /i is a smooth function. It 

must be kept in mind, however, that we cannot choose rj independently in the two regions 
r < ro and r > tq, because rj must be smooth at r = tq. As a consequence, we cannot 
conclude from ( |2.26| ) that and ^7 must both vanish on the Cauchy horizon. We only 
get the weaker condition that they have a similar behavior near this horizon; namely, it is 
necessary that the following "jump condition" holds: 

(7° + 7^) (^-A/lu+.ro + e,i?,vp+) - §/(n+, ro - e, i9, (^+) 
= o(l + |(7° + 73)§M(n+,ro + e,i9,(^+)|) as < e ^ 0. (2.27) 

On the event horizon, we obtain in the same way the condition 

(7° + 73) {^o{u+,ri+e,^,^+) - ^m(u+, n - e, (^+) 
= o(l + |(7° + 73)§Af(u+,ri-e,i9,v9+)|) as < e ^ 0. (2.28) 

The constructions we just carried out in the chart (U4., r, -i?, (/9+) can be repeated simi- 
larly in the coordinates (u_, r, ■(?, (^_). We list the resulting formulas: The transformation 
to the chart (m_, r, -i?, (/?_) gives for the Dirac matrices G^~ and G^~ 

To make these matrices regular, we transform the spinors according to ( p. 24 ) with 
V{r) = - (lAri + e(A) |A|3) 1 + 1 (\^\-\ - e(A) |A|3) 7V • 



2 

This gives for the transformed Dirac matrices 



a sin?? 2 ^"^ + 0^/0 •^n 

G"- = 



1 A)j^ - (1 + A)73 



2VU 
1 



& = — ^7' + ^(7° +7') 
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By evaluating the Dirac equation across the horizons in the weak sense, we obtain the 
conditions 



(7° - 7^) (§A/(w-,ro + - §/(n_,ro -e,i?,(^_)) 

= o(l + |(7°-73)§Af(n-,ro + e,i?,99_)|) as e - 

(7° -7^) (§o(u-,ri - §A/(u-,ri -e,i?,(^_)) 



0(1 + 1(7° - 7')^'Af(n-,ri - e,i},ip.)\) 



as £ 



(2.29) 



(2.30) 



It remains to transform the conditions (|2.27| )- (|2.3C| ) back to Boyer-Lindquist coordi- 
nates. Since the t- and (^-dependence of our wave function (2.17) has the form of a plane 
wave, we immediately conclude that the condition ( |2.27| ) is also valid in Boyer-Lindquist 
coordinates. We do not again write out that we consider the limit < e — > 0, but note 
that this condition was obtained by extending the Boyer-Lindquist space-time across the 
point (r = ro,t = —00), 

(7° + 73) (§M(t,ro + e,i?,v9) - ^ i{t,ro - e,!}, 



across t 



-00. 



(7° + 7^) ^((|Ar^ + |A|t) + 
|A| 



lAI 



lA 



^)7V) 



We next substitute the transformation ( 2.24 ) and ( |2.25| ) of the spinors. Using the identity 

1 
2 

(7° +7') 

we obtain the condition 

(7° + 73) {^M{t,ro+e,^,<f) 



|A| 



^/(t,ro -e,t?,v3)) 



0(1 + lAj-i 1(70 + q,^j(^t,ro + e,^,^)\) 



across t 



-00. 



(2.31) 



Finally, we consider the transformation ( |2.15| ). It also preserves the factor (7^^ -|- 7^) 
because 

(7° + 7')S 



diag y{r + ia cos t?) 2 , [r + ia cos t?) 2 , (r — ia cos ??) 2 , [r — ia cos t?) 2 

X |A|3 (7O -^7^) . (2.32) 

The diagonal matrix in this equation is irrelevant because it is regular on the horizons. 
The factors |A|^4 in (2.31) are compensated by the factor |A|4 in (2.32), and we end up 
with the simple condition 

(7°+7=') (l'Af(t,ro+e,i?,99) - ^'/(t, ro - e, (/;)) 

= o{l + \{-f^ +-f^)i!M{t,rQ + e,'d,ip)\) across i = -00. 

Similarly, the relations ( |2.28D -( 2.30D transform into 

(7°+73) (l'o(t,ri 99) - l'A/(t,ri -e,t?,(^)) 



0(1 + 1(7° + 7') ^Af(t,ri-e,i?,t)|) 



across t = 00 



(70-73)(^'A/(t,ro + e,i?,(/^) - ^'/(t, ro - e, i?, 92)) 



o(l + l(7°-7')^A/(i,^o + e,^?,¥')|) 



across t = 00 



(7O - 7^) {i>o{t, n + e, V>) - ^A/(i, n - e, i?, v)) 

= 0(1 + 1(7° -73)l'A/(t,ri -£,??, 99)1) across i 

The equations (|]3|)-(|!3|) are our matching conditions. 



-00. 



(2.33) 

(2.34) 
(2.35) 
(2.36) 
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Figure 1: The conformal structure of the non-extreme Kerr-Newman solution 
2.3 Non-Existence of Time-Periodic Solutions 

Before giving the proof of Theorem we must specify our assumptions on space-time 
and on the Dirac wave function in mathematical terms. By patching together the Kerr 
coordinate charts, one obtains the maximally extended Kerr-Newman space-time Q. We 
do not need the details of the maximal extension here; it suffices to discuss the Penrose di- 
agram of Figure |^. Abstractly speaking, the time and axial symmetry of the Kerr-Newman 
solution means that the metric admits two Killing fields. In order to better visualize these 
symmetries, one can isometrically map the regions of type /, M, and O into the regions 
r < vq, ro < r < ri, and r > ri of the Boyer-Lindquist coordinates, respectively; then 
the Killing fields are simply the vector fields dt and dip. The mappings to Boyer-Lindquist 
coordinates are unique up to the isometries of the Boyer-Lindquist space-time (i.e. ro- 
tations around the symmetry axis, time translations, parity transformations, and, in the 
case a = 0, time reversals). The Kerr coordinates, on the other hand, allow us to describe 
three adjacent regions of type /, M, and O including the boundaries between them. When 
we speak of Boyer-Lindquist or Kerr coordinates in the following, we implicitly mean that 
adjacent regions of the maximal extension are isometrically mapped to these coordinate 
charts. The maximal extension may be too general for a truly physical situation (if one 
thinks of a black hole which evolved from a gravitational collapse in the universe, for 
example, one may want to consider only one asymptotically flat region). Therefore we 
assume that our physical space-time is a given subset of the maximal extension. We call 
each region of type O which belongs to the physical space-time an asymptotic end. We 
assume that, in each asymptotic end, a time direction is given. Thus we can say that each 
asymptotic end is connected to two regions of type M, one in the future and one in the 
past. We assume that the physical wave function can be extended to the maximal Kerr- 
Newman space-time; in the regions which do not belong to the physical space-time, this 
extension \I' shall be identically zero. We want to consider a black hole, i.e. the situation 
where particles can disappear into the event horizon, but where no matter can emerge 
from the interior of a horizon. Therefore we assume that ^ is set identically zero in all 
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regions of type M which are connected to the asymptotic ends in the past. The remaining 
assumptions can be stated most easily in Boyer-Lindquist coordinates. Since the phase of 
a wave function ^ has no physical significance, when we say that ^ is time-periodic with 
period T we mean that there is a real parameter U such that^ 



^{t + T, r,^,ip) 



-iQT 



^{t,r,i3,<f) 



For time-periodic wave functions, we can separate out the time dependence in a discrete 
Fourier series. More precisely, we can write the wave function as a superposition of the 
form 



^{t,r,'d,(p) 



-int 



(2.37) 



with 



where the radial and angular functions X^^^{r) and Y^'^^{'d) satisfy the equations ( 2.20| ) 
and (|]2l|) with 

uj[n) = il + — n . 



The index A in ( |2.37D labels the eigenvalue of the operator A in (|2.19D ; the set cr^{A) 
denotes (for fixed n and k) all the possible values of A. As shown in the Appendix, the 
set <J^{A) is discrete. Finally, we specify our normalization condition: The Dirac wave 
functions are endowed with a positive scalar product (. | .). For this, one chooses a space- 
like hypersurface Ti together with a normal vector field u and considers for two wave 
functions ^ and $ the integral 



(^ I 



n 



^G^^ Uj dn 



(2.38) 



where ^' = ^'*7'^ is the adjoint spinor, and where djjL = ^ d^x is the invariant measure 
on 7i [g now denotes the determinant of the induced Riemannian metric). In a regular 
space-time, current conservation V j^G^^t = implies that the scalar product ( 2.3^ ) is 
independent of the choice of the hypersurface. The integrand of (^|^')>^ has the physical 
interpretation as the probability density of the Dirac particle. Therefore one usually 
normalizes the wave functions in such a way that (^' | ^) = 1. In our context, the 
singularities of the metric make the situation more difficult. Indeed, the normalization 
integral for a hypersurface crossing a horizon is problematic since such a hypersurface 
will necessarily fail to be space-like on a set of positive measure. We will therefore only 
consider the normalization integral in each asymptotic end away from the event horizon. 
More precisely, we choose for given r2 > ri the one-parameter family of hypersurfaces 



n 



t2 



{(t, r, ■(?, ip) with t = t2, r > r2} 



For a normalized solution ^ of the Dirac equation, the integral (vf | ^)'Ht gives the 
probability of the particle to be at time t2 in the region outside the ball of the radius r2 

^We remark that the condition of time-periodicity inside the horizon can be weakened to local uniform 
boundedness in t. This is proved by an "averaging argument" identical to the one given in jl^, Appendix A]. 
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around the origin; this probabihty must clearly be smaller than one. Therefore we impose 
for a normalizable solution the condition that, in each asymptotic end, 



I ^)nt, < oo 



for all t2- 



(2.39) 



We now begin the non-existence proof by analyzing the wave function in each asymp- 
totic end in Boyer-Lindquist coordinates. The following positivity argument shows that 
each component of the Dirac wave function must be normalizable: We average the nor- 
malization condition (|2.39|) one period and use the infinite series ( 2.37 ) to obtain 



oo 



1 /■*+^ 



n,n' k,k' A, A' 
X 



-i(k' —k)ip 



$Anfc(^^^) $Anfe . 



(2.40) 



Since the plane waves in this formula are integrated over a whole period, we only get a 
contribution if k = k' and n = n'. As is shown in the Appendix, the i?-integration gives 
zero unless A = A', and thus the right hand-side of ( p. 401) reduces to 



E E 

n,fc6Z AGo-^K-A) 



(2.41) 



Since the scalar product in ( 2.41| ) is positive, we conclude that the normalization integral 
must be finite for each ^^^^ 



[9 I <P < OO 



for ah n,k€Z, X€ a^{A), t G IR. 



In this way, we have reduced our problem to the analysis of static solutions of the Dirac 
equation. This simplification is especially useful because it enables us to work with the 
matching conditions of the previous section. In the following, we again use the notation 

Lemma 2.1 The function |Xp has finite boundary values on the event horizon. If it is 
zero at r = ri, then X vanishes identically for r > ri. 



Proof: The radial Dirac equation ( ^.20| ) gives for r > ri 



dr 



<i 



|A| —X, X> + <X, 
dr 



I A I d 
|A| —X> 

dr 



2ARe(X+X") + 2mrlm{X+X-) 



and thus 



(|A|+mr)|X|2 < ./|A| — |Xp < (|A|+mr)|X|2 
* dr 



(2.42) 



(2.43) 



In the case that |Xp(r) has a zero for r > ri, the uniqueness theorem for solutions of 
ODEs yields that X vanishes identically. In the opposite case |Xp(r) > for all r > ri, 
we divide ( 2.43| ) by ^| A| |Xp and integrate to get 



(|A|+ms)|A| 2 ds < log|Xp < / (|A| +ms) |A|~2 



ri < r < r . 
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Using that the singularity of |A| 2 at r = ri is integrable, we conclude that log |Xp has 
a finite limit at r = ri. ■ 

Combining this lemma with our normalization and matching conditions, we can show now 
that the Dirac wave function is identically zero outside the event horizon: 

Lemma 2.2 vanishes identically in each asymptotic end. 

Proof: According to our black hole assumption, we can apply the matching condition 
( 2.36| ) with = 0. Expressed in the radial functions X, this implies that 



lim X (r) = 

ri <r^ri 

Using asymptotic flatness and the transformation 



(2.44) 



the normalization condition (2.39) is equivalent to the integral condition 



\Xr dr < 00 



for r2 > ri. 



(2.45) 



As an immediate consequence of the radial Dirac equation ( 2.20| ), 



d_ 

dr 



X 



+ |2 



\x- 







(2.46) 



Thus the function - j^-p is a constant; the normalization condition ( |2.46D implies 

that this constant must be zero, 



= 



Together with the condition (|2.44| ) , we obta^iii that lini^-j^<^7- ^^--^ \x\ 

that X, and consequently niust vanish identically. 



0. Lemma 2.1 yields 



We remark that equation ( ^.461) can be interpreted physically as the conservation of the 
Dirac current in radial direction (notice that ^'C'^' = U~2 (|X+p — |X~p) |yp). 

It remains to show that the wave function also vanishes in the interior of the horizons. 



For this, we use the matching conditions and estimates similar to those in Lemma 2.1 



Proof of Theorem 1.1: According to Lemma 2.2, ^' is identically zero in all regions of 
type O. We first consider a region of type M in Boyer-Lindquist coordinates. Crossing 
the event horizon in the past and in the future brings us to regions of type O where ^ 
vanishes. Thus the matching conditions ( |2.34| ) and ( p.36| ) yield that 



lim (7° + 73)^'Af(t,r,??,^) = 

ri>r— >ri 



lim (7° -7^) ^'A/(i,r,??,^) 

ri>r— >ri 



and thus 



lim ^M{t,r,^,ip) 

ri>r^ri 



(2.47) 



The radial Dirac equation ( 2.2C| ) in the region M implies that 



|A|-1|X|2 
dr 
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(this again corresponds to the conservation of the radial Dirac flux). Thus |Xp(r) can 
only go to zero for ri > r — > ri if it is identically zero. We conclude that ^ must be 
identically zero in all regions of type M. 

Finally, we consider a region of type I in Boyer-Lindquist coordinates. We can cross 
the Cauchy horizon at t = cxd or t = — oo; this brings us to regions of type M, where 
= 0. Thus the matching conditions ( |2.33| ) and ( |2.35 ) imply that 



lim ^ lit^r^-d^ip) 

ro<r^ro 







(2.48) 



The radial Dirac equation ( p.20D inside the Cauchy horizon is the same as in the asymp- 
totic end. Thus the estimate ( |2.43| ) again holds, and we conclude from (2.48) that 
vanishes identically. ■ 



3 The General Separable Case 



Our goal in the present section is to prove an analogue of Theorem LI for the most general 
stationary axisymmetric and orthogonally transitive metric in which the Dirac equation 
can be separated into ordinary differential equations by Chandrasekhar's procedure^ An 
expression for this metric was determined in jlj]. When the Einstein-Maxwell equations 
are imposed, this metric gives rise to all the generalizations of the Kerr-Newman solution 
discovered by Carter as well as to a family of exact solutions for which the orbits of 
the two-parameter Abelian isometry group are null surfaces |^]. 

Our aim will be to state as many of the hypotheses of our theorem as possible in purely 
geometric terms. We will thus begin by giving a geometric characterization of the metric 
which constitutes the starting point of We will limit ourselves to the case when the 
orbits of the isometry group are time-like 2-surfaces, since the procedure is quite similar in 
the case of space-like or null orbits. We wish to point out that the assumption of orthogonal 
transitivity is a very natural one to make in a general relativistic context. Indeed, it is a 
classical theorem of Carter and Papapetrou ^ that if the energy-momentum tensor 
satisfies some mild invariance conditions, then every stationary axisymmetric solution of 
the Einstein equations has the property of admitting through every point a 2-surface 
which is orthogonal to the orbit of the isometry group through that point, i.e. the metric 
is orthogonally transitive. 

Let us first recall from that any four-dimensional metric of Lorentzian signature 
admitting a two-parameter Abelian group of isometrics acting orthogonally transitively 
on time-like orbits admits local coordinates {u,v,w,x) in which 



ds' = T 



-2 



{Ldu + Mdvf - {Ndu + Pdvf - ^ - ^ 



(3.1) 



where the metric coefficients L, M, N, P and the conformal factor T ^ are functions of 
X and w only, and where X = X{x), W = W{w). Furthermore, we have LP — MN ^ 0, 

■^We should point out that there exist metrics admitting only one Killing vector, in which the Dirac 
equation is completely separable into ordinary differential equations jll|. These metrics are not covered 
by Chandrasekhar's approach since the symmetry operators underlying this different type of separability 
are necessarily of order 2 or higher. In contrast, Chandrasekhar's method always gives rise to symmetry 
operators of order 1 since it fits into Miller's theory of factorizable separable systems jtL5|. 
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T > 0, > 0, X > in order for the metric to have the required Lorentzian signature. 
It is manifest that the action on the isometry group generated by the Kihing vectors 
^ and ^ is orthogonally transitive since the orbits, which are given by the time- like 
2-surfaces x = ci, w = C2, admit the orthogonal 2-surfaces given hy u = c'l, v = 
Carter [Q] proved furthermore that the four-dimensional Lorentzian metrics admitting 
an Abelian isometry group acting orthogonally transitively on non-null orbits have the 
following remarkable property: there exists in the isotropy subgroup at every point an 
element of order two whose differential C maps every vector X in the tangent space to the 
orbit through that point to its opposite, and maps every vector Y in the normal space to 
the orbit to itself. In the local coordinates {ujVjW, x), this involutive isometry is given by 
mapping ,it;,x) — > (— n, — f Recall now that a Newman-Penrose null frame for 
the metric is said to be symmetric ^ if under the involution £, we have 

£(/) = —n , ; C{n) = —I , £(m) = — m , £(m) = —m . (3-2) 

A symmetric null frame {l,n,m,fn) for the metric (3.1) is given by 



m 



LP 


- MN 




1 


LP 


- MN 




1 


LP 


- MN 




1 


LP 


- MN 



ou ov ) aw 

ou ov J ow 



m = —i — ^ ( -M — + L — - tVX — 

ou ov I ox 



T f 1 / d d\ ^ d 

ou ov I ox 



Besides the metric ( p.lD , we will also consider the electromagnetic vector potential given 
by 

A = Ajdx^ =T{U{Ldu ^ Mdv)^ K{N du ^ Pdv)) , (3.3) 

where = H{w) and K = K{x). The 1-form (|3.3| ) is manifestly invariant under the action 
of the isometry group of the metric ( |3.lD . Note that the Maxwell field 2- form F = dA 
is skew-invariant under the involution C, so that the electromagnetic energy-momentum 
tensor is invariant under in accordance with the hypotheses of the Carter-Papapetrou 
theorem (01. 



It was shown in [14| that a necessary condition for the Dirac equation (G — m) ^ = 
to be separable in the above symmetric null frame for the metric (^]^) and in the 
electromagnetic vector potential (|3.3D is that the metric functions L, M, N and P satisfy 
the constraints 

d L ^ _ 9 M 

dx^LP - Mn' ~ ' dx^LP - MN' 
d N _ d P 

dw^LP - Mn' ~ ' dw^LP - MN' 

In geometric terms, these conformally invariant conditions are equivalent to the require- 
ment that the flows of the null vectors fields / and n of the symmetric frame be geodesic 
and shear-free, 

PVjl^ = , 2{Vjk + ^kh) - {V%)g,k = , 
n^Vj-n^ = , 2(Vjnfc + V^nj) - (V^'np)^^^ = . 
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It is noteworthy that these conditions are necessary and sufficient for the Stackel separa- 
bihty in the metric of the Hamilton-Jacobi equation for the nuh geodesic flow, 

ikdS dS 
9 = 

By implementing these separability conditions and doing some relabeling of the metric 
functions and the coordinates we see (see ||9|, Section 3] ) that the metric can be put in the 
form 



W Z X z 

— (ei du + m dv)"^ — — dw^ — — {e2 du + p dv)"^ — — dx^ 
Z W Z J\ 



, (3.4) 



where ei and 62 are constants and where 

m = Tn{x) , p = p{w) , Z = eip — . 

The next necessary condition for the Dirac equation to be solvable by separation of vari- 
ables in the symmetric frame for the metric ( |3.lD and the electromagnetic potential ( |3.3| ) 
is is that the algebraic structure of the Weyl conformal curvature tensor be of type D in 
the Petrov-Penrose classification [14|. This conformally invariant condition is necessary 



and sufficient for the separability of the massless Dirac or Weyl neutrino equation = 



see [14, Theorem 2]). The type D condition for the metric ( |3.4D can be expressed by the 



condition that the 1-form u given by 

UJ := (AZ)^^ {eim' (x) dw + e2p'{w)dx) 

be closed, 

du = 0. (3.5) 

We therefore have locally 

00 = dB 

for some real- valued function B{w,x). 

Finally, it was shown in |l^. Theorem 3] that the massive Dirac equation (G — m)"^ = 
is separable in the symmetric frame for the metric ( |3.4| ) if in addition to the conditions 
stated above, there exist real-valued functions g{x) and h(w) such that 

^1/2 

T ex.p{2iB) = h{w)+ig{x) . (3.6) 

This is the only separability condition for the Dirac equation which is not conformally 
invariant. 

We begin by analyzing the type D condition |3.5| . It was proved in [^, Section 4] that, if 
this condition is satisfied, then a set of coordinates {u, v, w, x) for the metric |3.4| can be 
chosen in such a way that the metric functions p{w),m{x), the constants ei,e2 and the 
function B take one of the following four forms: 

Case A: 

ei = 1 = 62 , p{w) = vj^ , m{x) = —x^ , 

i (w-ix 

4 \w + IX 



18 



Case B_ 



ei 
B 







62 



k 



4 \K + tx 



p{w) = 2A;it; , m{x) 



Case i?4 



ei 



1 , 62 = 

i f w — il 
4 



m(x) 



-2lx 



Case COO; 



ei 



1 , 62 

-7r/4 . 



p{w) 



m(x 



(x) = , 



We are using the same labelHng for the cases given above as the one used in ||l|, |^. Each 
of these cases admits an invariant characterization which we will not recall here but which 
is given in |jl^. Next, by imposing the separability condition ( |3.6| ), we conclude that the 
conformal factor must be constant in each of the cases A, B^, and C^^. This 
constant can be normalized to one by a rescaling of the coordinates. 

We will restrict our attention for the rest of this section to Case A, which can be thought 
of as the generic case. Indeed, we will see below that the solutions one obtains when 
imposing the Einstein-Maxwell equations in Case A contain the Kerr-Newman solution as 
a special case. On the other hand, it is shown in |]|, that the solutions one obtains when 
imposing the Einstein-Maxwell equations in the remaining cases have isometry groups of 
dimension 4 in Cases B^, Bj^ and 6 in Case C^^ . These solutions notably include the Taub- 
NUT and Robinson-Bertotti metrics. Furthermore, it is shown in the above references that 
all these solutions can be obtained from Carter's A solution by a suitable limiting process. 

The metric in Case A is thus given by 



W{w) 
w"^ + x"^ 



{du — x'^dvY 



VJ^ +x^ , 9 



W{w) 



X{x) 

^2 _|_ J.2 



{du + uP' dvY 



up' + 



dx'^ . (3.7) 



The Dirac equation for the metric (|3.7|) can be separated by a procedure analogous to 



the one described in Section 2 for the Kerr-Newman metric. It is proved in [14| that the 
spinor field given by ^ 



^{u, V, w, x) 



/ {w -ixf/"^ X-{w)Y-{x) \ 
{w - ixfl'^ X^iw) Y^{x) 
{w^ixfl^ X^{w) Y-(x) 

\ {w ^ ix) 



(3.8) 



^ To facilitate the comparison with the expression given in (2.15) for the Kerr-Newman case, we have 
chosen to include the factor 14^3 in the transformation of the spinor field. Just as in the Kerr-Newman 
case, this transformation has the effect of producing a slightly simpler eigenvalue equation in the variable 
vj than the one obtained in lll4|l. 
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where a; and k are constants, will be a solution of the Dirac equation expressed in the 
symmetric frame for the metric ( |3.7| ) and the Weyl representation of the Dirac matrices, 
if and only if the transformed spinor 



^(w, x) 



( x^{w) y_(x) \ 

X+{w)Y+{x) 

x+{w) y_(x) 

V X^{w)Y+{x) ) 



satisfies the eigenvalue problems given by 



where 



X 



\ 











im w 













—im w 


^ 









—im w 





V 







im w 


—m X 





^ 


fx Cx+ 





m X 


-VxCx_ 





y 


/X Cx-[. 


—m X 





—VX £x~ 








m X 



(3.9) 



(3.10) 



and Piu-i- and Cx± are the ordinary differential operators defined by 

1 



d_ 

dw 
d 



W{w) 
1 



dx ^ 'X{x) 



-iujw + ik — ieH{w) 



iujx^ — ik — -X'{x) — ieK{x) 



As a motivation for our generalization of Theorem 1.1, we recall from how the Kerr- 
Newman solution arises as a special case of the metric (3/7) and vector potential (|3.3| ). 
For this, we first impose the Einstein-Maxwell equations 



Ri 



1 



1 



-Rgij - A gij = FikF^ - -QijFuF 



kl 



(3.11) 



which determine the remaining functions X, W , and the electromagnetic field J- = dA. 
The general solution is given by 



W 
X 
A 



A 
3" 

A 4 
3" 
1 



W'' + f2W^ + fiW + /o + + 



f2X +gix + fo 



(^Q w {du — dv) + P X {du + vo^ dv)j 



(3.12) 
(3.13) 
(3.14) 



where Q and P denote the electric and magnetic monopole moments and /o,/i,/2i5i 
are arbitrary parameters. These solutions were first discovered by Carter The Kerr- 
Newman solution is obtained as a special case by putting suitable restrictions on the pa- 
rameters A, /o, /i, /2) ffi appearing in the metric and by implementing through the choice 
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of coordinates (u, v, w, x) the additional hypothesis that the metric is stationary axisym- 
metric (as opposed to simply admitting a two-dimensional Abelian isometry group). This 
procedure will serve us as a guide in formulating the conditions under which the analogue 



of Theorem 1.1 holds for the Dirac equation in the metric (|3.7| ). 

We first re-label the coordinates (u, v, w, x) as {t, ip, r, ii) and rescale the ignorable 
coordinates t and (p so as to have ei = 1, 62 = a, where a is a constant. Note that while a 
is freely normalizable at this stage, it will become an essential non-normalizable parameter 
once our freedom to scale the coordinates will have been exhausted. The expressions of 
the metric functions p{r) and m(/x) will thus be slightly different in this normalization of 
the coordinates. We have 

p{r) = , m(^) = — a/i^ . 

Taking the cosmological constant A to be zero, the metric functions W'^{r) and X'^{fi) 
appearing in Carter's A solution will take the form 

W{r) = f^r^-2Mr + foa^ + Q^ + P^ , X{fi) = -f2 fx^ + gi fx + fo ■ 

We now implement the hypothesis of axisymmetry of the metric by letting ip be an 
angular polar coordinate adapted to the axis of symmetry of the metric and taking /x to 
be proportional to the cosine of the polar angle measured from the axis. The range of (p 
will thus be the interval (0, 27r) and the range of ^ will be a bounded open interval. In 
order for the metric to be of hyperbolic signature, we must require that X{^) be positive 
for varying in that bounded interval. It follows that the roots of must be distinct 

and that X{^) must be positive as 11 varies in the interval bounded by these roots. This 
gives the constraints 

/2 > , gl + 4f2fo > . 

Next, we can ensure that the singular behavior of the metric at the roots of X{fi) is 
caused by nothing more than the usual angular coordinate singularity at the polar axis 
by choosing the ignorable coordinates and the remaining parameters appearing in the 
metric in such a way that the roots of m(^) and X(fj,) coincide. Thus, we use our residual 
freedom to replace t by a constant linear combination of t and cp to add the same arbitrary 
constant to p{r) and m{^). The roots of m{^) will then be located at the endpoints of a 
symmetric interval (— c, c). In order for the roots of X{^) to also be located at c and — c, 
the coefficients gi , /2 and /o must necessarily satisfy the relations 

51 = , /2 > , /o > . 

We can now set /2 = /o = 1 by rescaling ^ and r, so that the range of becomes the 
interval (—1, 1). With these normalizations, the parameter a has now become an essential 
parameter in the metric. The metric functions p, m, W and X thus reduce to 

p{r) = r"^ + , m{fj,) = a{l — , 

W{r) = -2Mr + a^ + , X{^i) = 1 - /U^ . 

We see that by letting cos ^ = 1), where — vr < -i? < tt, we recover the Kerr-Newman metric 
in Boyer-Lindquist coordinates (|2.l| ). 

We are now ready to state a theorem which generalizes Theorem 1.1 to the most 
general family of stationary axisymmetric metrics in which the Dirac equation is solvable 
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by separation of variables. We will consider a normal form for the metrics (|3.7| ) of Case A 
in which the stationary and axisymmetric character of the metric is made manifest through 
an appropriate choice of coordinates and by the imposition of suitable restrictions on the 
singularities of the metric functions. This normal form is easily established by a procedure 
similar to the one we described above for the Kerr-Newman metric. 

Theorem 3.1 Consider the stationary, axisymmetric metric and the vector potential 
given by 

= ^2 + ^2^2 {dt -a{l-fi) dip) ^^^^ dr 

^^^^ {a dt - (r2 + a2) d^)^ - + d/.^ , (3.15) 



A 



(Hir) (dt - a (1 - ^2) dip) + K(fi) (a dt - (r^ + a^) dip)) , (3.16) 



where a > 0, -(x, < t < (x,, < p < 2tt, H{r) G C°°(IR), and K{p) G C°°([-l,l]). 
Assume that the functions X{fi) £ C°°([— 1, 1]) and W{r) £ C°°(IR) have simple zeros 
at = 1, —1 and r = ri,. . . ^r^, so that the range of the coordinates is fi £ (—1, 1) and 
r G IR\ {ri, . . . , rTv}- Suppose furthermore that the metric is asymptotically Minkowskian, 

< lim W{r) < oo . 

r— »itoo 

Then the Dirac equation (G — m) \I' = has no normalizable time-periodic solutions. 
This theorem can be thought of as an analogue in the axisymmetric context of the gen- 



eralization presented in |13, Remark 4.3] of the non-existence theorem for normalizable 
time-periodic solutions of the Dirac equation in the Reissner-Nordstrom background. We 
can likewise argue here that the proof of Theorem is similar to that of Theore m |1.1| . 
Indeed, we can see from (|3.1C| ) that the separated equations obtained for the metric (|3.15| ) 
are similar in structure to those obtained in the Kerr-Newman case. Next, we remark that 
in view of the assumptions made in Theorem 3.1 to ensure the regularity of the metric at 
the symmetry axis, it follows that the maximal analytic extension of the metric ( |3.15| ) will 
have a conformal diagram similar in structure to the one obtained for the non-extreme 
Kerr-Newman metric (see Figure ||). We thus conclude that matching conditions for the 
spinor fields across the horizons take a form identical to the equations ( |2.33| )-( ^.36| ) which 
were obtained for the Kerr-Newman case. The key observation is then that the assumption 
that we made about the zeros of W{r) will guarantee that the estimate required to prove 
the analogue of Lemma 2.1 is valid in this more general context. Finally, the regularity of 
the eigenfunctions of the angular equation is established by the same procedure as in the 
Kerr-Newman case, by showing that the angular operator can be viewed as an essentially 
self-adjoint elliptic operator on the 2-sphere S'^ with C°° coefficients. 

We conclude by remarking that non-existence theorems similar to Theorems 1.1 and 3.1 
likewise hold true in Cases 5+, S_, and C^^. 



A Regularity of the Angular Part 

In this appendix, we will show that all solutions of the angular equation ( |2.21| ) are regular. 
More precisely, we will see that these functions are of class in the open interval 
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< "i? < TT and uniformly bounded on the closed interval < ■!? < vr. The method is 
to reduce the problem to an elliptic eigenvalue equation on the 2-sphere, where standard 
elliptic regularity theory can be applied. 
Consider the PDE 



ia -— H + ia — — + acv sin-i^u — am cos?? a ] a = Xa (A.l) 

ydtf 2 J smv dip J 

for a two-component, complex function a{'&, ip),0<'&<'ir,0<^p< 2tt, with the boundary 
conditions 

lima{'d,ip) = — lim a{'d,ip) 

ip\0 



{p^ are again the Pauli matrices (2^)). With the ansatz 



a 



the eigenvalue equation (A.l) simplifies to ( 2.21| ). Thus it suffices to show regularity for 



the solutions of ( A.l 



Unfortunately, the coefficients in ( [A.lD have singularities (which are not just removable 
by a coordinate transformation on S"^). However, after performing the transformation 



a — > a = \Ja with lJi(d^Lp) = exp^— i — a J exp^— i — cj 
we obtain for a the equation Aa = Ad with the smooth operator 
A = i \ cjV — o"' ) + aoj sin aip — am cos ■!? 



dr 

where a'^ and a"^ denote the Pauli matrices in polar coordinates, 

= sin"!? cos(/9 0"^ + sint? sin 990"^ + cos'd 

LP ■ 1 I 2 

= — sm ip a + cos ip a 

Note that the transformation U is not continuous on 5"^, as it changes sign along the 
line 99 = 0. Combined with our above boundary conditions for a, we conclude that 
d is a continuous function on 5^. Hence we can regard the equation Aa = Ad as an 
elliptic eigenvalue equation on the compact domain S"^. The operator A is essentially self- 
adjoint on C°°(S'^)^ C L^(5^)^. Thus its square is a positive, essentially self-adjoint 
operator with smooth coefficients. Standard elliptic theory yields that has a purely 
discrete spectrum with finite-dimensional eigenspaces and smooth eigenfunctions. Since 
the eigenvectors of A are obtained by diagonalizing A on the finite-dimensional eigenspaces 
of A^, they are also smooth. 
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